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Abstract 

Background: Peste des petits ruminants (PPR) is an endemic and highly contagious disease in small ruminants of 
Pakistan. Despite the fact that an effective vaccine is available, outbreaks are regularly occurring in the country. Thus 
so far, the diagnosis has primarily been made based on clinical outcome or serology. This study was carried out to 
characterize PPRVfrom an emerging wave of outbreaks from Punjab, Pakistan. 

Results: A total of 32 blood samples from five different flocks were tested with real-time PCR for the presence of 
PPRV genome. The samples detected positive in real-time PCR (n = 17) were subjected to conventional PCR for the 
amplification of the nucleoprotein (N) gene. Phylogenetic analysis of the sequenced N genes (n = 8) indicated the 
grouping of all the sequences in lineage IV along with PPRV strains from Asian and Middle East. However, 
interestingly sequences were divided into two groups. One group of viruses (n = 7) clustered with previously 
characterized Pakistani isolates whereas one strain of PPRV was distinct and clustered with Saudi Arabian and 
Iranian strains of PPRV. 

Conclusions: Results demonstrated in this study expanded the information on the genetic nature of different PPRV 
population circulating in small ruminants. Such information is essential to understand genetic nature of PPRV strains 
throughout the country. Proper understanding of these viruses will help to devise control strategies in PPRV 
endemic countries such as Pakistan. 



Background 

Peste des petits ruminants (PPR) is a highly contagious viral 
disease of domestic and wild ruminants where it cause high 
morbidity (100%) and mortality (90%) [1,2]. The causative 
agent, peste des petis ruminants virus (PPRV), is grouped in 
the genus Morbillivirus within family Paramyxoviridae along 
with rinderpest, measles, phocine-, dolphin-, canine- and 
porpoise-distemper viruses [3]. Like other Morbilliviruses, 
PPRV is pleomorphic, with negative sense single stranded 
RNA genome containing 15,948 nucleotides [2,4]. The gen- 
ome follows the "rule-of-six" and encodes six structural pro- 
teins that include nucleocapsid (N), phosphoprotein (P), 
matrix (M), fusion (F), haemagglutinin (H) and large poly- 
merase (L) [4]. 

Based on the molecular characterization, strains of 
PPRV can be grouped into four lineages, which are 
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genetically distinct from each other. Lineage I includes 
isolates from Western Africa, lineage II contains isolates 
from West African countries, the Ivory Coast, Guinea 
and Burkina Faso, and lineage III represents strains from 
Eastern Africa, the Sudan, Yemen and Oman [2,4-6]. 
The lineage IV comprises of PPRV strains from the 
Arabian Peninsula, the Middle East and South Asia 
[2,4,6]. Classification of PPRV is being analyzed based 
on the sequence analysis of both F and N genes; how- 
ever, parallel comparison of PPRV strains has proposed 
that N gene is most divergent and hence most appro- 
priate for molecular characterization of closely related 
isolates [7]. The virus has been recognized to occur as 
only one serotype among four lineages [6]. 

PPRV has been documented in Pakistan since 1991, 
however, the isolates were confirmed through PCR few 
years latter in 1994 [8]. Since then, PPR remained en- 
demic in Pakistan despite use of a live attenuated vac- 
cine in small ruminants. Due to serological monitoring 
facilities, it remained difficult to ascertain the level of 
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vaccine failure and thus aggravates disease epidemiology 
and its control Therefore, determining the nature of cir- 
culating strains in different parts of Pakistan is crucial to 
not only aid in disease diagnosis and but also to devise 
better control strategies in future. The present work has 
been conducted to determine the genetic nature of cir- 
culating PPRV strains that are continuously causing out- 
breaks in central Punjab, Pakistan. 

Methods 

The study involves one of the livestock richest districts of 
Punjab where an emerging wave of clinical disease, 
suspected for PPRV, was reported. A brief clinical history 
and outcomes of different diagnostic tests are demon- 
strated in Table 1. The number of animals in the herd and 
their ages varied from 27-50 animals and 3 months- 
5 years, respectively. The breeds of animals were either 
Beetal or non-descript, without any history of vaccination 
in the past. Food and water were provided ad libitum. 
The morbidity and case fatality rate was in the range of 
50-80% and 20-35%, respectively. A majority of infected 
animal died within first week of infection. Whole blood 
samples (n = 32) were collected aseptically from animals 
from five different outbreaks. The blood (200-300 ul) was 
poured onto QIAcard FTA Indicator Four Spots (Qiagen, 
Hilden, Germany), which preserve genomic material and 
lysed the cells and viruses. 

The total RNA was eluted from the QIAcard FTA Indi- 
cator as described [9] and was used to screen by real-time 
PCR for the presence of PPRV genome as reported earlier 
[9,10]. All the samples appear positive in real-time PCR (Ct 
values <35) were used in conventional PCR to amplify N 
gene of PPRV for downstream sequencing [9,11]. The re- 
sultant PCR products were gel extracted and processed for 
sequencing using ABI PRISM BigDye Terminator version 
3.1 (Applied Biosystems), according to the manufacturers 



instructions. Sequences were analyzed with an automated 
nucleic acid analyzer (ABI PRISM 3100; Applied Biosystems). 
Each DNA fragment was sequenced at least twice in both 
directions. 

The sequences were edited and assembled using EditSeq 
and SeqMan suits within Lasergene 8 (version 8.0.2 13, 
DNASTAR, Inc., Madison, WI, USA). The resultant se- 
quences were aligned with the sequences retrieved from 
GenBank representing all the lineages of PPRV. Construction 
of phylogenetic trees was performed with the neighbour- 
joining method using Kimura two-parametermodel in Mega5 
version 5 (CEMI, Tempe, AZ, USA). 

Results 

Out of total samples (n = 32) analyzed, real-time PCR 
detected only 18 samples with Ct values lower than 35. 
The results of real-time PCR confirmed the cause of 
these outbreaks to be PPRV. However, to demonstrate 
the genetic nature of PPRV, a conventional PCR for the 
amplification of N gene was performed. Although all the 
samples were not detected positive, a representative 
sample of each outbreak was amplified which appear 
sufficient for sequencing and downstream analysis. All 
the N genes were submitted to the GenBank under ac- 
cession number KC207867-KC207874. 

Nucleotide sequence identity among eight sequences 
reported here and with previously characterized Pakistani 
strains of PPRV was found to be 97.2-100% and 94.6-97%, 
respectively. Phylogenetic tree, based on the 255 bp of the 
N genes sequenced in this study and collected from 
GenBank, was shown in Figure 1. All the Pakistani strains 
of PPRV clustered into lineage IV, which is exclusively 
prevalent in Asian and Middle East countries. However, it 
was interesting to note that all the Pakistani strains of 
PPRV, sequenced here or reported earlier, were grouped 
closed to Tajikistani strains of PPRV, except one sequence 



Table 1 Brief history of outbreaks and outcome of different diagnostic assays 


Place of sampling 
in district Okara 


Clinical findings 


Animals sampled/ 
total animals in 
the herd 


Real-time PCR 
(positive/total) 


Conventional PCR 
for N gene 
(positive/total) 


Chak 48/2-L 


High fever (105-107 F), nasal discharges, sneezing, 
coughing, sticking of mucus to nostrils, erosions in the 
oral mucosa. Severe diarrhea and then drop in temperature 
after start of the diarrhea. Death nearly 5-10 days after 
start of disease 


8/35 


3/8 


1/8 


Chak 54/2-L 


Fever (104-105 F), nasal discharges, sneezing, coughing, sticking 
of mucus to nostrils, erosive ulcerative stomatitis, diarrhea 


3/27 


3/3 


1/3 


Chak 1/44. 


Fever (104-105 F), nasal discharge, sneezing, coughing, necrotic 
stomatitis, severe diarrhea 


9/39 


5/9 


2/9 


Chak 31/24 


Body temperature (105-107 F), weakness, nasal and lacrimal 
discharges, sneezing, coughing, erosive lesions in the oral cavity, 
severe diarrhea 


4/50 


3/4 


1/4 


Chak 24/2-R 


Fever (104-105 F), purulent nasal discharges, coughing, erosive 
ulcerative stomatitis, diarrhea 


8/41 


3/8 


3/8 
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Figure 1 Phylogenetic analysis of the PPRV. Neighbor-joining majority rule consensus tree from variable N gene (255 bp) region constructed 
with Kimura two-parameter model in Mega5 version 5 at bootstrap value 1000 replicates. In the figure bootstrap values only above 50% are 
shown. The horizontal lines were proportional to the distance among sequences and the samples characterized in this study were marked with 
black-square (■). 



(Pakistan/Okara/MM61/2012). This strain was clustered 
with PPRV isolates from Saudi Arabia and Iran, and col- 
lectively made a separate cluster within lineage IV. The 
Paldstan/Okara/MM6 1/2012 strain carried eight amino 
acid differences from rest of the strains studied here. 
Taken together, the results demonstrated that there are at 
least two different population of PPRV circulating in the 
country. 

Discussion 

Since first report of PPRV in Pakistan, several outbreaks 
have been reported [2,9,12-19]. Most of these were based 
on either clinical history, antigen or antibody detection. 
However, molecular characterization of the virus is 



considered a key to understand the genetics of the 
causative agents. The results of real-time and conven- 
tional PCR from blood samples, transported onto FTA 
cards, confirmed the fact that they are suitable, reliable 
and economical mean of transportation [9,19], especially 
where maintenance of cold chain is not possible or legis- 
lations for shipment of biological material is complicated 
[2]. Conventional and real-time PCR for the detection of 
PPRV have shown variable sensitivity and specificity. 
Therefore these tests must be evaluated in parallel be- 
fore their implementations for routine diagnostic pur- 
poses in specific laboratory settings. 

Latest data has shown that N gene reveals improved 
image of PPRV epidemiology and is considered better 
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than F gene based distribution of PPRV strains [6,20]. As 
this fact is well established that PPRV has propensity to 
change genetically, N gene based phylogenetic analysis 
was performed in this study. The similarity and diver- 
gence among samples from five outbreaks and formerly 
characterized PPRV isolates were evaluated by analyzing 
their common N gene sequences. In general, the top- 
ology of the tree revealed PPRV isolates from diverse 
geographical areas differ largely in their partial N gene 
sequence. The tree shows that currently analyzed PPRV 
are more closely related to Tajikistan. The possible rela- 
tionship of Pakistani strains to that of Tajikistan strains 
of PPRV could be due to the fact that Pakistan with its 
province named Gilgit Baltistan shares the Pamir region 
with Tajikistan through Wakhan corridor with Afghanistan 
where animal movement in the mountainous grazing lands 
of Pamir region is not uncommon between these countries 
(Figure 2). Moreover, selling of small ruminants by nomads 
in Gilgit Baltistan province to rest of the country, on or be- 
fore religious ceremonies such as Eid-ul-Adha, is common 
particularly in Punjab province that represents about 60% 
of total population. The grouping of these strains of PPRV 
at two different places in the tree shows that PPRV isolates 



may have undergone exceptional evolutionary pathway or 
may be a result of several introduction of viruses from a 
variety of sources. However, such interpretations required 
confirmation by the analysis of larger samples collected 
from wide range of outbreaks. 

Conclusions 

In conclusion, lineage IV of PPRV is currently circulat- 
ing in the country, with certain level of genetic diversity. 
Since serological assays (i.e. cELISA) don't necessary in- 
dicate the current persistence of the infection, it is es- 
sential to implement molecular diagnosis along with 
characterization. Despite PPRV is endemic in Pakistan 
and outbreaks are regularly occurring, limited informa- 
tion is available on the genetic nature of PPRV. The se- 
quences reported here expand the available information 
on the circulating strains of PPRV in Pakistan. There is 
no sequence information for any of the N genes of PPRV 
is available other than our previously characterized se- 
quences and the sequences reported here. With these 
limited data, we still manage to demonstrate the pres- 
ence of two different population of PPRV, which warrant 
future need to perform such studies at country level to 




Figure 2 Detection of PPRV in Pakistan and surrounding countries, a) Pakistan; b) Iran; c) Afghanistan; d) Turkmenistan; e) Uzbekistan; 
f) Tajikistan; g) China; h) India; i) Nepal. Dates shown on the map illustrate both published reports and reports to the OIE reference laboratories as 
detailed in Banyard et al., 2010 (20). For Pakistan, provinces are shown: /') Gilgit baltistan/Northern Areas; /'/') Khyber pakhtunkhwa; /'/'/') Tribal areas; /V) 
Punjab; v) Azad Jammu & Kashmir; vi) Baluchistan; w7) Sindh province. The approximate position of the outbreak detailed in this study is 
represented by a star. + India frequently reports outbreaks as PPRV is endemic across the country (20). 
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ascertain the complete picture of circulating viruses in the 
country. Such understating is crucial for devising future 
control plans. Currently, vaccination is recommended in 
certain areas of the country. This vaccination is based on 
Nig75/1, which belong to lineage II, while field isolates 
from Pakistan are grouped in lineage IV. Genetic 
characterization of field strains will provide founda- 
tions for construction of vaccines from domestic 
strains as has recently been practiced in India. 
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